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1. Introduction 



The meander problem is one of tiiose tantalizing questions that has resisted a definite 
solution for decades, although it is very easy to state: "Given a positive integer n, in how 
many topologically distinct ways can a closed non-intersecting planar curve {road) cross a 
straight line [river) in exactly 2n points {hridges)T^ 

Originally an exercise of recreational mathematics the meander problem turned 
out to have applications in the most various branches of science: Sorting algorithms in 
computer science 0, enumeration of ovals of planar algebraic curves P], classification of 
three-manifolds |^ , and in connection with a particular type of self-avoiding walk describing 
the compact folding of a linear chain ^j. 

An obvious strategy would of course be to evaluate the first few meander numbers 
Mn, in the hope of finding an explicit formula, valid for arbitrary n. Such enumerative 
approaches exist on various levels of sophistication ||^-|^, and a recent transfer matrix 
method |jTy] carried out this program up to n = 24. Although an explicit expression for 
Mji appears to be out of reach, it became clear that, in analogy with two-dimensional 
lattice polymers, the meander numbers scale asymptotically as M„ ~ CR'^^/n°', where R 
is a connectivity constant and a a configuration exponent. 

A major achievement of random matrix theory has been to deal with precisely such 
asymptotic enumeration problems . It is therefore natural to apply such techniques to 
the meander problem HT^ - IT^ . In particular it has emerged that a generalized multi-road 
multi-river meander problem, in which each closed segment of river (resp. road) is given 
the statistical weight ni (resp. 712), can be cast as a Hermitian matrix model, known as the 
0(^1,712) model 0. In the special case of 712 = 1 this model is soluble by a saddle-point 
method, leading to an exact evaluation of R and a [[T5| , p!6[] . Unfortunately these results do 
not pertain to the original meander problem, which is recovered in the limit ni, 712 —>■ 0. 

In a recent publication [|T^ it was argued that the meander problem is a particular 
realization of the coupling to gravity of a certain two-fiavored loop model [jl^, initially 
defined on the square lattice. The most general gravitational version of this loop model is 
a generalization of the meander problem in which river and road segments, counted with 
their respective weights of ni and 77-2, are allowed to cross as well as to touch one another 
without crossing (tangency points) [[l^ . We shall refer to this model as tangent meanders. 
On the regular lattice, directed segments of river and road can be inserted by means of 
certain magnetic defect operators, for which the anomalous dimensions are known exactly. 
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When dressed by quantum gravity, these dimensions transform according to the KPZ 
formula |T^ . This transformation aUows one to extract exact values for the configuration 
exponent a of tangent meanders, whereas R, being a non-universal quantity, is lost in the 
process. The connection to the original meander problem is then made by arguing that 
tangency is irrelevant from a renormalization point of view |T^. Thus, the result for a in 
fact pertains to the original meander problem, i.e. to the gravitational 0(ni,n2) model. 
Moreover, the operator content of the theory gives access to other geometries, involving 
several rivers, possibly with marked points, as well as semi-meanders (river with a source). 

Here we review and extend the arguments of [^. In particular we establish the 
irrelevance of tangency rigorously in a number of special cases. We also add credibility to 
the theoretical predictions by performing extensive exact enumerations of various meander 
geometries, using a generalization of the transfer matrix method presented in [|10| . We first 
address multi-component meanders, that allow for precisely checking the predicted value of 
the central charge of the underlying conformal theory. Next we explore two distinct river 
geometries, namely (i) two parallel rivers, and (ii) one semi-infinite river, that permit to 
validate their magnetic operator formulation within the corresponding conformal theory. 
We finally consider the case of tangent meanders and verify the irrelevance of tangency, 
thus confirming the cornerstone of the argument. In all cases we find an excellent agreement 
with theory, typically confirming the configuration exponents with 4-5 significant digits. 

The paper is organized as follows. In Section 2 we review the square-lattice loop model 
and its solution, before presenting its gravitationally dressed version and the results for 
the asymptotics of a range of meander-related quantities. The transfer matrix algorithms 
are presented in Section 3, and in Section 4 we analyze the data and compare them to the 
theoretical predictions. Our conclusions and some perspectives can be found in Section 5. 



2. Theory: from fully-packed loop gases to meanders 



In this Section, we review the arguments of Ref. |17] relating the meander problem 
to the gravitational version of a particular fully-packed loop model initially defined on the 
square lattice |]T^ . The effect of gravity is to replace the lattice with a random quadrangu- 
lation of the sphere. The lattice loop gas is described in subsection 2.1 while its conformal 
structure is presented in subsection 2.2. The connection to meanders via two-dimensional 
quantum gravity is explained in subsection 2.3. This leads naturally to an effective field 
theory description detailed in subsection 2.4 together with the subsequent predictions for 
various meandric configuration exponents. 
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Fig. 1: A typical fully-paclced loop configuration on the square lattice. 
Assuming doubly periodic boundary conditions, there are 6 black loops (solid 
lines) and 4 white ones (dashed lines). Up to rotations, the vertices of the 
model are of the two types (a) "crossing" or (b) "avoiding" . 

2.1. Fully-packed loop gas on the square lattice 

The configurations of the fully-packed loop model that we shall consider are defined 
by assigning to each edge of the two-dimensional square lattice either of two colors (say, 
black or white, represented as solid or dashed lines in Fig. 1, also referred to as 1 and 2 in 
the following) in such a way that each vertex has exactly two black and two white incident 
edges. Up to obvious rotations, this gives rise to the two vertex configurations depicted 
in Fig. 1 in which the black and white lines either avoid or cross each other. Note that 
with periodic boundary conditions the black and white lines form loops. It is interesting 
to remark that the fully-packed loop model's configurations defined here differ from those 
of the so-called densely packed loop model [jl^ in that each vertex is visited by a black 
and a white loop, whereas in the dense case, loops of a given (say black) color are not 
constrained to visit all vertices. 

The partition function of the fully-packed loop model is then defined by assigning a 
weight ni per black loop and n2 per white one, 



^FPL(ni, n2) 



E 



(2.1) 



fully — packed loop 
configurations 



where we have denoted by Li the total numbers of loops of each color i = 1, 2. Following [r^| 
we shall denote the model with partition function ( |2.1| ) as the FPL^(ni, 77-2) model, while 
the densely packed version is referred to as the DPL^(ni,n2) model. 

The loop weights may be recast as local Boltzmann weights as follows. This step 
is important in the field theoretic description of the model, since it leads to a local field 
theory. Let us assign to each black or white loop an arbitrary orientation, and attach to 
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each vertex a local Boltzmann weight e™^^^^^^'^'^^'^)!'^ where = 1 if the oriented loop of 
color i makes a left turn, = if it goes straight, and = — 1 if it makes a right turn. 
Summing over all possible orientations of all loops, we get a factor 2 cos ttc^ per loop of 
color z, and therefore we reproduce the desired loop weights by setting 



ni = 2 cos TTCi 



77.2 = 2 COS 7re2 . 



(2.2) 



2.2. Conformal field theory description 




-> © A 



■< — O B 



o C 
O D 



Fig. 2: A typical configuration of the FPL^ model together with the bicol- 
oration of its vertices (checkerboard of filled (•) and empty (o) dots). We 
have added the corresponding dictionary that allows to map the loop config- 
urations onto A, S, C, D labelings of the edges. 



The FPL^(?7,i, model is known to be critical for < < 2 [jT^, and is described in 
the continuum limit by a simple conformal field theory based on free scalar fields. To iden- 
tify its basic degrees of freedom, it is useful to rephrase the model as a (three-dimensional) 
height model as follows. Starting from an oriented fully-packed black and white loop con- 
figuration, we first bicolor the vertices of the square lattice, say with alternating filled (•) 
and empty (o) dots. Then we use the dictionary of Fig. 2 to assign one of the four labels 
A, B, C, D to each colored and oriented edge. With this convention, it is clear that edges 
of type ABAB . . . alternate along black loops, whereas edges of type CDCD . . . alternate 
along white loops, and that each vertex has one incident edge of each type A, B, C and 
D. It is seen that the four-labeling with A,B,C,D is in one-to-one correspondence with 
the coloring and orientation of edges of the FPL'^ model. In particular, the orientation of 
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Fig. 3: Rules determining the change of the height variable across la- 
beled edges. We adopt the Ampere convention that the height is increased 
(resp. decreased) by the edge value if the arrow of the edge points to the 
left (resp. right). The edge labels must be interpreted as three-dimensional 
vectors with the respective values A, — B, C, — D. 

a given black or white loop is reversed if we interchange the A and B or C and D labels 
along the loop. 

The above colors allow for the definition of a dual vector height variable on the center 
of each face of the lattice. Indeed, viewing as vectors the A, B, C, D labeling of the edges 
of the lattice, let us arbitrarily fix the height to be zero on a given face of the lattice, 
and define it on all other faces by successive use the rules of Fig. 3 for the transition 
from a face to any of its neighbors. Note that it is necessary to impose the condition 
A-l-B + C + D = 0to ensure that the heights are consistently defined around each vertex. 
We may therefore assume in all generality that A, B, C, D are actually four vectors in IR^ 
with vanishing sum. To get a more symmetric formulation, we may further fix A, B, C, D 
to be the four unit vectors pointing from the center of a tetrahedron towards its vertices. 
The heights are then clearly three-dimensional, as linear combinations of A,B,C,D. In 



the continuum limit, it was argued |T^ that the three-dimensional height variable turns 
into a three-dimensional scalar field. Moreover the symmetries of the model completely 
fix the action for these fields and the corresponding field theory is conformal, with central 
charge 

e2 



CFPL(ni,n2) = 3-6 -^ + -^ , (2.3) 

V 1 - ei 1 - 69 / 



ei 1-62. 

where Ci have been defined in ( p.2|) and are constrained by < 6^ < 1/2. Actually the 
shift in the central charge away from 3 is due to the introduction of a background electric 
charge, ensuring that loops that have non-trivial winding with respect to the periodic 
boundary conditions still get correctly weighted, although for such loops the argument 
given before (|2.2|) no longer holds true. 



2.3. Meanders: the coupling to gravity 

To finally get to meanders, we must consider the coupling of the FPL^(ni, 77.2) model to 
two-dimensional quantum gravity, by allowing the square lattice to fiuctuate into arbitrary 
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planar four-valent graphs. For each such graph the fuUy-packed loop model is still defined 
by coloring the edges black or white and allowing only the vertices shown in Fig. 1. As 
before, each colored loop is weighted by the appropriate factor (z = 1,2). 

If we try to go through the steps of the previous section, namely by transforming the 
model into a height model, the issue of bicolorability of the vertices of the lattice becomes 
crucial on a random four-valent graph. Indeed, not all such graphs are vertex-bicolorable. 
So the coupling to gravity stricto sensu (sum over arbitrary planar four-valent graphs) will 
destroy this property. 

We may now follow either of the two following paths. First, we can repair this and 
impose that the particular coupling to gravity preserve the bicolorability, namely that the 
gravitational model be defined on the set of vertex-bicolorable four-valent graphs only. 
These graphs are dual to the so-called Eulerian quadrangulations. In genus zero (planar 
case), the latter are characterized by the fact that all their vertices have an even valency 
(Euler condition). If we couple the FPL^(ni, ^2) model to Eulerian gravity, the A, B, C, D 
labeling is still well-defined, as well as the three-dimensional height, now defined on the 
centers of the faces of the graph. This preserves the degrees of freedom of the fiat space 
model entirely in the gravitational formulation. This approach was initiated in [EDI for the 
simpler case of the fully packed 0(n) model with only one type of loops. 
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Fig. 4: Rules determining the change of the height variable across labeled 
edges in the non-bicolored case. These rules are identical to those of Fig. 3, 
with the further restriction that B = —A and D = — C, allowing to ignore 
the bicoloration of vertices. 



On the other hand, we can study how the degrees of freedom of the model are affected 
by the coupling to ordinary (non-Eulerian) quantum gravity. Having lost the bicolorability 
of vertices, it is no longer possible to distinguish between A and B labels on one hand, 
and C and D on the other. We may still define an edge-labeling of the graph in one-to-one 
correspondence with oriented colored fully-packed loop configurations on the graph, but 
with vectors A, B, C, D satisfying the two constraints A + B = and C + D = 0, and 
picking A and C to be two perpendicular unit vectors in IR^ . The correspondence between 
color /orientation and A, C labels reads as in Fig. 4. As in Eq. ( |2.1|) , the model is further 
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completed by attaching weights to each loop of color i = 1,2. We may now define a 
height variable on the centers of the faces of the graph, by use of the previous rules. The 
main difference is that the height now lives in two dimensions (the plane generated by A 
and C). Such a dimensional reduction is also observed on the square lattice when going 
from the FPL^(ni,n2) to the DPL^(ni, 77.2) model ||18|| ; interestingly enough, the same 



dimensional reduction is also observed when the FPL^ model is defined on the Manhattan 
square lattice, with oriented loops respecting the Manhattan orientation [^|. It results in 



a shift c ^ c — 1 in the central charge of the underlying conformal theory, namely 



c(ni,7i2)=CDPLK,n2)=2-6 -^ + -^ . (2.4) 

VI -ei 1-62/ 

In the following, we will concentrate on this formulation, eventually leading to the solution 
of the meander problem. 

The partition function of the fully-packed model coupled to ordinary quantum gravity, 
hereafter referred to as the GFPL^(ni,n2) model, reads in genus zero: 

^GFPL(ni,n2;x,y)= Yl E |Aut(r c)! ^^'^^^^'^-^^^?/'^^^'^^ . (2-5) 

four — valcnt planar FPL configs.C ' \ ' / I 

graphs r on T 

where the sum extends over all the genus zero four-valent graphs F, and |Aut(F, C)| is 
the order of the symmetry group of F equipped with the loop configuration C. We have 
also denoted by 14, Vb the total numbers of vertices of type a and b defined in Fig. 1 in 
the particular loop configuration, namely we have weighted each crossing of a black and 
a white loop by x and each avoiding by y. When x = y, these are interpreted as the 
cosmological constant, as the total number of vertices Va + V;, of F is also the area of the 
corresponding dual random surface. 
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Fig. 5: Examples of (a) a tangent meander with 2 tangency points (b- 
vertices) and 6 bridges (a-vertices) and (b) a meander with 8 bridges. 
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To generate meanders, we must now extract from (2^) the configurations with only 
one black and one white loop, that will respectively play the role of the river and the road. 
This is done by taking the limit ni, ^2 ^ in (|2.5|), resulting in 



ZGFPL{x,y)= lim (ZGFPL(ni, ^2; y) - 1) 

ni.na^O nin2 



^-^ x^^yP (2.6) 

n,p>0 ^ 

7T,-|-p> 1 



where we have denoted by /i2n,p the total number of tangent meanders with In crossings 
and p tangency points, i.e. configurations of a non-selfintersecting circuit (road) crossing a 
line (river) through 2n points (bridges) and touching the river p times (tangent contacts), 
as illustrated in Fig. 5. The usual meander numbers defined in correspond to only 
crossings and no tangent points and read = fX2n,o- In (|2.6|), the prefactor l/(2(2n + p)) 
stands for the symmetry factor attached to the tangent meanders: going from the graph 
to the representation where the river is a line, we may indeed cut the river loop in 2n + p 
places in between bridges and tangency points, and we still have 2 choices for the up/down 
position. 

NE NE 
h+A| h h+AV h 



h fh-C ^ h .A-C 

SW SW 

Fig. 6: A type b vertex of the FPL^ gravitational model, together with its 
dual height configuration. We note that the NE and SW heights are identical. 
We may therefore undo the vertex as shown, which results in its irrelevance. 

The meanders are therefore generated by the function (|2.6| ) for y = 0. Let us now 
show that the universality class of the tangent meanders is the same as that of meanders. 
In the transformation into a (two-dimensional) height model, the "tangency" vertex b of 
Fig. 1 corresponds to the arrangements of heights on adjacent faces depicted in Fig. 6. We 
notice that the NE and SW heights are identical, irrespectively of the orientations of the 
two loops. This means that as far as the height variable is concerned this vertex may be 



simply removed as shown. We conclude that the b- vertex of the model is irrelevant [17|. As 



a consequence, we expect the following asymptotics for N = 2n + p large of the partition 
function //Ar(x,y) for tangent meanders with a total of bridges and tangency points: 

l,N{x,y)= Yl x'V/^2n,p~C(x,y):^^||^, (2.7) 

2n+p=N 

where the configuration exponent a is independent of x and y for x > and y > 0. In 
particular, when y = 0, a is identified with the meander configuration exponent M„ ~ 
cK^/n"'. The irrelevance of the vertex b will be checked numerically in subsection 4.4 
below. 

More generally, we may consider tangent meanders with one single river, but arbi- 
trarily many possibly interlocking closed roads. This is given by the tt-i ^ limit of ( |2.5| ) 
while n2 = q is kept finite, namely 

ZGFPL{q;x,y) = lim — (ZcfplI^i, ^2 = x, y) - 1) 

"-1^0 m 

rj,2rLyP (2-8) 
n,p>0;n+p>l ^ ^' 

where we have defined the tangent meander polynomial 

n+p 



^2n,p(?) = 5^^£V (2-9) 

fc=l 



with coefficients iJi'2np being the numbers of tangent meanders with k connected com- 
ponents of road, 2n bridges and p tangency points. The polynomial /i2n,o(?) = "^n(Q') 
coincides with the meander polynomial defined in f^. We may also define the canonical 
partition function: 

fiN{q;x,y)= J2 a;2"yP^2n,p(?) • (2.10) 

Tl ,p > 

2n+p=N 

We expect a large asymptotic behavior of the form 

m(g; X, y) ^ C{q- X, y) ^''^^(^^^ , (2.11) 

where the configuration exponent a[q) only depends on q and not on y/x, due to the 
irrelevance of the b- vertex. In particular, it takes the same value at y/x = 0, where 
it coincides with the multi-component meander configuration exponent, i.e. mn{q) ~ 
C(g)i?(Q)2'^/n"(9). 
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In the special case q = 1 of arbitrarily many roads without extra weight, these num- 
bers can be computed exactly. Indeed, we may decompose an arbitrary multi-component 
tangent meander with 2n bridges and p tangency points into its upper part (above the 
river) and lower part (below) , and consider that these may be obtained first by picking pi 
points among the total = 2n -|- p to be the upper tangency points, and p2 = p — pi to 
be the lower ones. Let us then draw small semi-circles tangent to the p points, pi of them 
in the upper half, p2 in the lower. With the 2n bridges, we have now a total of 2n + 2pi 
points in the upper half and 2n + 2p2 in the lower one to be connected among themselves 
by pairs through non-intersecting arches. Such upper and lower arch configurations have 
already been extensively studied in 0. In particular, there are Cm = (2m)!/(m!(m+ 1)!) 
[cm are the Catalan numbers) distinct upper arch configurations connecting 2m bridges 
by pairs in the upper half plane above the river. Hence we must choose among the c^+pi 
upper arch configurations and the Cn+p2 lower ones to form an arbitrary multi-component 
tangent meander. This results in the formula 

PiKp - PiV.{^n)l 

where the combinatorial factor accounts for the choices of upper and lower tangent points 



among the total of 2n + p. The corresponding partition function (|2.10|) reads 

2;2n ypi yP2 

fiN{q = l;x,y) = Nl ^ ^ ——cn^^^cr,^^^. (2.13) 

2n+p=N Pl+P2=P 

When is large, this is easily estimated by a saddle-point technique, making use of the 
Stirling formula. For y > 0, we find the large behavior 

^^(Q = l;x,y)~^^^±^ (2.14) 

up to a multiplicative constant depending on x and y only. This shows explicitly that the 
exponent cx{q = 1) = 3 is robust and is not affected by the respective values of x and y. 
This confirms in particular the above-mentioned irrelevance of the b- vertex, as a keeps to 
the same value, irrespectively of y. 

In conclusion, the meander numbers belong to the universality class of the GFPL^(?t,i, 77-2) 
model at 711,^2 0. The corresponding fiat space theory has the central charge ( p.4|) with. 



according to ( |2.2|) , ei = 62 = 1/2, hence c = —4. In the next subsection, we will concen- 
trate on the partition function Zgfpl(2^, 2/ = 0) ( |2.6| ) that generates the meander numbers. 
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More generally, the multi-component meander polynomial belongs to the universality 
class of the GFPL^(ni, 712) model, with ni ^ and 122 = q and central charge ( p.4|) . 
Finally, we may also consider meanders with arbitrarily many rivers and roads, with a 
weight ni per river and 77-2 per road; these objects belong to the universality class of 
the GFPL^(ni,n2) model. It has been shown using matrix model techniques that the 
GFPL^(?7.i = 1,^2 = q) with y = belongs to the same universality class as the 0(g) 
model when coupled to ordinary gravity (i.e. defined on arbitrary four-valent graphs) 
|15|j . More precisely, the number of multi-river, multi-road meanders with a total of 2n 
intersections and with a weight q per road and 1 per river, also weighted by their inverse 
symmetry factor (and multiplied by 4n to make it comparable to Mn) behaves for large n 
as 

C(l, ,)2"/n"(i.»l, R{l,q) = 2 ""''7^^' , 

a{l,q) = , g = 2cos(7re). 

1 — e 

The determination of R{l,q) on which we have no prediction was made possible by an 
explicit mapping of the matrix model onto a particular version of the gravitational 0{q) 
model, where the critical value of the cosmological constant Xc can be explicitly calculated. 
The value a(l, q) corresponds to the expected scaling behavior of a c = 1 — 6e^/(l — e) 
(c.f. (|2.4| ) with TT-i = 1 and 122 = q) conformal theory coupled to two-dimensional quantum 
gravity. In particular, when g = 1, c = we get a(l, 1) = 5/2 as expected in "pure gravity" 
without matter. 



y — r ^ \-- y 




Fig. 7: The transformation of the four-valent a- and b-vertices of the FPL^ 
model into pairs of trivalent ones connected by an extra (gray) edge, depicted 
as a wavy line. Due to the 1 to 2 correspondence in the case of the a-vertex, 
we must restrict ourselves to the = 2y, and then each trivalent vertex 

receives a weight ^/y. 
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It is interesting to notice that if we restore the b- vertex (with a weight y per vertex), 
and choose specificaUy x = 2y, then the GFPL^(l,(j) model can be mapped onto that 
of tricoloring the edges of a random trivalent graph (with black, gray and white edge 
colors) with a weight q per loop of alternating say gray and white colors. Indeed, the 
quartic vertices a and b may be decomposed into pairs of connected trivalent ones as 
shown in Fig. 7 in which we have added a third type of edge, say with gray color. This 
transformation allows to map the configurations of the GFPL^(l,q') model weighted by 
X = 2y per a-vertex and y per b-vertex onto those of arbitrary edge-tricolored trivalent 
graphs, with a weight ^/y per vertex and q per white/gray loop. The latter model was 



solved in and identified with a particular version of the fully-packed 0(2^) model on 
random trivalent graphs: indeed, the white and gray edges form naturally fully packed 
loops on the graph, and on each loop we may interchange the two colors to generate new 
tricolorings, hence we may simply draw fully-packed loops on trivalent graphs, and attach 
a weight 2q to each of these loops. But the loops have all even lengths and this turns out 
to dimensionally reduce the model to an effective 0(2^/2 = q) one with central charge 
c = 1 — 6e^/(l — e) as before. This particular point therefore lies in the same universality 
class as the model without b-vertex, which gives another explicit example of the irrelevance 
of the b-vertex. 

2.4- Field theory description of meandric numbers 

The coupling of a conformal field theory with central charge c < 1 to two-dimensional 
quantum gravity (i.e. its definition on random surfaces) has a simple field-theoretical 
formulation in terms of the Liouville field describing the conformal classes of metrics of 
the surfaces. This has led to a number of results, including the precise determination of 
various critical exponents. Indeed, the gravitational theory (say on genus zero surfaces) 
displays a critical behavior as a function of the cosmological constant x. In particular, 
there exists a finite value Xc of x at which the (connected) partition function behaves as 

Z(x) ~ (xc-x)2-^(^), (2.16) 



where the string susceptibility exponent 7 is related to the central charge c through [|T9 



.(0 ^ -'-v-';-"^^-' (2.17) 
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When applied to the DPL2(0,0) model of the previous section (with c = —4), whose 
gravitational version was shown to describe meanders, we find that 



7 = 7(-4) = ^ . (2.18) 



Comparing (|2.16| ) with the expansion 



^2n 



^gfpl(x) = J2 -^^r. (2.19) 

n>l 



we deduce the asymptotic behavior [17 



M„~C^, „ = 2-,= H^^. (2.20) 

n° ' ' 12 ^ ' 

Moreover, a number of the operators of the flat space conformal theory (in particular 
the spinless ones, with conformal dimensions h = h) get dressed by gravity, in such a way 
that they acquire anomalous scaling dimensions. Any given operator (pk with dimensions 
hk = hk, is dressed into an operator cpk with dressed dimension Afc such that the correlation 
functions behave as 

. ..4>k,) ~ (x, - x)2-^+^(^'=.-i) (2.21) 

when X approaches the critical value Xc, and where the dressed dimension is related to 
the flat space conformal dimension hk through 

Let us now present the operator content of the c = —4 conformal theory describing 



the dense loop model DPL (0, 0) [T^]. For generic values of ni,n2, the DPL (77,1,77-2) has 
a continuum description as a two-component scalar field with charges at infinity. More 
precisely, it is a Coulomb gas made of two decoupled scalar fields, with c = c{ni) = 
1 — 6ef/(l — ei) = — 2 at = (e^ = 1/2) respectively, each viewed as the effective field 
theory of loops of one color. In particular, within each scalar field theory (indexed by the 
color z = 1, 2), there exist operators '(l)^j^\z) that create k oriented defect lines (of color i) 
for the scalar field, with conformal dimensions ||23|| : 



(,) ^ ft 4 

fc 32 V ; 
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at Hi = 0. In the Coulomb gas formalism, these correspond to electromagnetic opera- 
tors with electric charge (spin- wave) and magnetic charge ±k/2 (vortex), according to 
whether the defect line is oriented from or to the insertion point, and /c is a strictly posi- 
tive integer. The electric charge ensures that, if the defect lines wind around the insertion 
point, all extra curvature weights get cancelled For /c = 0, (|2.23| ) must be replaced 
by h^^ = corresponding to the identity operator. The correlation functions must have 
a vanishing total magnetic charge. Although these operators also carry an electric charge 
1/2 = 61=62, the electric neutrality of correlators imposes no extra condition0. 

We may also combine operators for both colors, namely consider mixed operators 
'4'ki,k2 = "^ith conformal dimension /ia;i,A:2 = h^k} + h^k2 ^i? ^2 G 2- 

Let us now study the dressing of these operators by gravity, and interpret them in 
meandric terms. The dressed operator -0^ again corresponds to the creation of a vertex 
with \k\ outcoming {k > 0) or incoming (/c < 0) lines of color i. 

As a first example, note that = 0. From the formula ( |2.21| ), it is clear that the 
operators ip^2 have the effect of marking a point on a loop of color i. However, due to 
the constraint of global magnetic neutrality of correlators, the marking is not arbitrary 
as the loop of color i must have an even number of marks (alternating 1^2^ and "ijj^^l)- In 
particular, the one-point function {ijj^^) that would naively count meanders with a marked 
river turns out to vanish, while for instance the two-point function {^jj^^ indeed counts 
meanders with two marked points on the river. We will meet more examples of this below. 



Fig. 8: A typical semi-meander configuration. 

^ It can be shown [p5| that the corresponding non-minimal conformal theory pertaining to 
the loops of color i (with central charge c = 1 — 66^/(1 — Ci)) is equally well defined as a theory 
with central charge c = 1, in which the operator content has been reorganized, and in particular 
all of the above electromagnetic operators get a vanishing electric charge, while their magnetic 
one is unchanged. A supersymmetric version of this theory was also introduced in [ p6[ , in which 
arbitrary correlators could be calculated, irrespectively of the electric charge. 
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The first appfication of the above concerns the two point function describing the 
insertion of a segment of river (color 1) 

(2.24) 

RecaU that when we take the hmit ni,n2 — > 0, only diagrams with one connected compo- 
nent of river and one of road are selected. In the case of (|2.24|) , the river forms a segment, 
around which the road can freely wind (see Fig. 8). To fix ambiguities, let us send one end 
of the river to infinity (say to the left) and therefore represent the river as a half-line (this 
is allowed as we work on Riemann sphere). The number of configurations of a closed road 
crossing a half-line (river with a source) through n bridges is defined as the semi-meander 
number M^. We immediately identify the series expansion of ( p.24| ) as a function of x to 
be 



(2.25) 



n>l 



We therefore deduce the semi-meander asymptotics 

~ a = l + 2Ai -7 = 1 + 



24 



(75 + ^29). 



(2.26) 




(a) 




■♦ * ^— • — f ^ ♦■ 

> • ^ <- 



(c) 

Fig. 9: A typical example of an eight figure river geometry (a). On the 
sphere, it is equivalent to the situation of (b), where the river crossing has 
been cleared of all winding pieces of road. This crossing may finally be sent to 
infinity (c) so as to form two parallel rivers. We have represented a particular 
configuration of road with 2pi = 4 bridges on one loop and 2p2 = 2 on the 
other. Due to magnetic neutrality in the operator formulation, the two loops 
are marked as shown. 
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Our second example is the "eight" meander, in which the river forms the figure "eight" . 
To generate it, we need to insert an operator '04^'*, but magnetic neutrahty forces us to 
insert two river-marking operators 'ijj^J^l, one on each loop of the eight. The three-point 
correlation function ('*/'4^'*('0-2)^) therefore generates the numbers of meanders whose river 
forms an eight figure, and with one marked point on each loop (see Fig. 9 (a)). To fix 
ambiguities, let us send the river-crossing to infinity. We may now represent the river as 
two parallel lines (all connected to the four-valent point at infinity), but the markings are 
still there (c.f. Fig. 9 (c)). Hence we have 



n>l 

4e 



^2n^2-mark 



X 



2n 



n>l 



J2 (2pi + l)(2p2 + l)Mp, 

Pi+P2=n 



(2.27) 



P2 5 



where "^"^^ is the number of meanders with two rivers, each of which is marked, and 



Pl,P2 



is the number of meanders with two parallel rivers, and with 2pi bridges one the 
first and 2p2 on the second (there are indeed 2pi + 1 possible markings on a river loop with 
2pi bridges) . Note the prefactor of 1 /2 accounting for the fact that we distinguish the top 
and bottom of the figure when we represent a two-river meander. From (|2.21|) , we find 



that 17 



2 — mark 



const. - 



X, 



-2n 



Q2-mark 



n 



a2-mark = A4 + 2A2 - 7 = ^(V^+ V^)(V5 + V^)- 



(2.28) 



Fig. 10: A meander with two parallel rivers may be viewed as a meander 
"with a seam". The point where the seam is attached is a remnant of the 
point at infinity (say to the left), where the two rivers meet. The seam just 
prevents roads from encircling this point. 



16 



It is also easy to determine the large n behavior of the unmarked meanders with two 
parallel rivers. Let us draw the two rivers horizontally. We may view the two rivers as 
connected to each other at their far left (at infinity) and represent them as in Fig. 10 as 
one marked line, such that a seam originates from the marked point that prevents any road 
from encircling the point to the left. So we have mapped two- river meanders onto one- river 
meanders with a seam. We may now generate all meanders with a seam by considering 
all the ways of placing a seam on the upper-half of meanders. If E{M) denotes the total 
number of upper exterior arches of a meander M (namely arches that have no other arch 
above them), then there are E{M) -\- 1 distinct ways of decorating the meander with a 
seam. Hence the total number of two-river meanders reads 

Mn~""= {E{M) + l) = {E + l)^Mr^, (2.29) 

meanders M 
with 271 bridges 

where (. . .)„ stands for the average over all meanders with 2n bridges. 




Fig. 11: How to construct a new meander with 2n -|- 2 bridges from a 
meander with 2n bridges and one distinguished exterior arch, here on the 
upper half of the meander. We cut the corresponding arch, and glue it back 
across the river, by encircling the whole lower part of the meander. This is 
easily inverted within the set of meanders with 2n -|- 2 bridges and one unique 
upper or lower exterior arch, showing that this construction gives only rise to 
distinct meanders. 

We now argue that {E + 1)„ remains finite when n — > oo. Indeed, starting from a 
meander with 2n bridges, we may construct one new meander with 2n + 2 bridges for each 
upper and lower exterior arch. As shown in Fig. 11, we simply cut the arch and close it on 
the other half, thus creating two new bridges and keeping a single connected road. Such 
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meanders have 2n + 2 bridges and one single exterior arch on the side where we have closed 
the arch, and are all distinct. Hence we may write that their number is bounded by M^+i, 
namely 

meanders AI ^ 
with 2n bridges 

We therefore have the asymptotic bounds 1 < (-E + 1)^ < 1 + -R/2, which implies that 
M^~^^^- oc R'^'^/n'^ like ordinary meanders, hence a2-riv. = «• We also deduce that 

2^Pi+P2=n^^-^Pi,P2 (2.31) 

13 = a- a2-mark = 2 - A4 = 2 - ^{VU -V5){V5 + V29) = 1.521898 • ■ ■ . 

This shows the rather unexpected result that the two-river meanders tend to be very 
asymmetric, with a number of bridges of the order n on one river (since pi + P2 = n) and 
of the order n^~^ on the other, with /3 — 1 = 0.521898 ■ ■ ■. The above arguments can be 
generalized to the case of multi-component meanders and give rise to similar predictions. 
In [T^, a number of other results have been presented, all corresponding to more 



sophisticated river geometries, and making use of the magnetic defect operators i^f. ■ It 
would also be possible in principle to use mixed operators to generate diagrams with both 
road and river geometries fixed. Another direction consists in going away from the point 
ni = 77-2 = 0, for instance by considering the GFPL'^(ni = 0, n2 = q) model in which 
meanders with arbitrary numbers of roads are considered, with a weight q per road. In 
that case, the corresponding conformal theory with central charge 



c=-l-6--^ — e = -Arccos(^) (2.32) 
1 — e TT 2 

may be viewed as two decoupled bosonic field theories: a c = —2 theory (that of river 
loops, at ni = 0) and one with c = c{q) = 1 — 6e^/ (1 — e), where q = 2 cos ire (that of road 
loops). The generating function Z'gfpl('?; = 0) ( |2.8| ) for meander polynomials reads 

ZGFPL(g; x,y = 0) = J2 -^'^r.iq) ~ (xc - a;)2-^('=) (2.33) 

n>l 

with 7(c) as in ( |2.17| ), and c as in ( p.32| ). This turns into the asymptotics 



n 



2n 



^^M-Ciq)-^, a(Q) = 2-7(c), (2.34) 
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namely [|T7 



1 - e + 3e^ + v/(l - e + 3e2)(13 - 13e + Se^) 
= 2 + — -^^ (2.35) 



with e as in (|2.32|) . An analogous formula holds for the (multi-component) semi-meander 
polynomials, namely 

~ ^(^)^^' ^(5) = 1 + 2Ai(c) - 7(c) (2.36) 



with Ai(c) as in (^^, c as in ( ^^3^ ), and hi = -3/32 as in (^]2|). This yields [|T7 



^ , V2(24e2 + e-l)(Vl-e + 3i^+Vl3-13e + 3e2) 
= 1 + ^^^^ (2.37) 

with e as in ( |2.32| ). 

An important remark is in order concerning the range of validity of ( p.35|) and ( |2.37|) . 
First, the DPL^(0,q') model is critical only for q < 2, i.e. e > 0. We expect therefore a 
very different scaling behavior for the meandric numbers when q > 2. At large q, it was 
shown in Ref. that 



a{q) = ^ (2.38) 



independently of q. 



Fig. 12: A typical "branched" semi- meander. 

On the other hand, for semi-meanders, this transition at g = 2 is never reached. 
Indeed, another phenomenon appears at a value qc < 2. It corresponds to a proliferation 
of connected components of road and also of the pieces of road that wind around the 
source of the river. More precisely, let us compare the two following contributions to the 
multi-component semi-meander numbers, according to whether: 
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(i) "branched" semi-meanders dominate: the semi-meander is typicaUy separated into 
two parts, namely the interior and exterior of a "big circle" as shown in Fig. 12. The 
contribution of such semi-meanders reads asymptotically 

^ mg^ Rjqr 

a(<j) ^ a(g) ~^2a(g)-l- ^^"^^^ 
ni+n2=n ""1 ""2 

(ii) "connected" semi-meanders dominate: the semi-meanders cannot be cut as in (i) , and 
typically behave as 

(2.40) 

^a(g) ^ ^ 

The transition between the two regimes (i)-(ii) will take place when 2a{q) — 1 = a{q). 
We deduce that precisely at the transition, we must have a = 1, which according to ( |2.37| ) 
takes place when 24e^ + e — 1 = 0, namely at the critical value e = q = qc given by 



97-1 „ V97-1 

2cos7r^^— , (2.41) 



48 ' " 48 

also corresponding to c = 3/4. Hence the formula (|2.37| ) is only valid for Cc < e, namely 
q < qc = 1.673849 • • • < 2. Beyond q = qc, we expect the average number of connected 
components of road to be of the order of n, and we have 

a{q) = for Q > q^. (2.42) 

It is quite interesting to notice that the transition of semi-meanders is of a different 
nature than that of meanders. The latter simply encounters the q = 2 transition of the 0{q) 
model. The former undergoes a winding transition, in which semi-meanders become very 
different from meanders. Indeed, as long as the dominant semi-meanders have very little 
winding around the source of the river, they behave just like meanders (hence R{q) = Riq) 
in the regime q < qc)- But when the winding number of the semi-meanders becomes 
relevant (of the order of n) we expect many more semi-meanders than meanders with the 
same number of bridges, and R{q) > R{q) for q > qc- 



3. Meander enumeration algorithms 

In this Section, we will describe the algorithms that we used to check numerically our 
predictions of Sect. 2 for the different configuration exponents. All the results presented 
below concern the case ni = 0, i.e. a connected river configuration, and a varying weight 
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n2 = q per road, i.e. multi-connected road configurations. We will also consider the several 
meander geometries discussed in Sect. 2. 

All our algorithms are based on transfer matrix techniques. Transfer matrices have 
proven to be very powerful for studying a wide range of statistical mechanics systems, 
especially in two dimensions. Originally applied to systems with local interaction, such 
as the Ising model, the "row to row" transfer matrix describes the transition between 
two successive rows of spins, say at 'time' t and t + 1 along the transfer direction. More 
recently, the transfer matrix technique was also applied to systems with non-local degrees 
of freedom, but where, in a suitably designed basis, the statistical weights can still be 
evaluated between neighboring time slices. In this way it became possible to study self- 
avoiding polygons and walks [^, and the random cluster model [^, to mention but 
a few important examples. 

Finally, transfer matrices can also be applied to random lattice problems, as long as a 
definite transfer direction can be defined, as for instance in the case of Lorentzian gravity 
||30|| . This is also precisely the case in the meander problem, where we can simply choose 
to transfer along the river, adding one bridge in each time step. This was first recognized 



by Jensen in where the method proved to be much better than previous enumeration 
algorithms, leading in particular to the largest accessible numbers of bridges. 

When implemented on a computer, the transfer matrix algorithms will allow to enu- 
merate exactly the various meandric objects that we are interested in for a fixed finite 
number of bridges up to a certain maximal value of N (typically up to 48 bridges in 
the results presented below). From these exact finite N values, we can extract estimates 
for the large N asymptotic behaviors, and in particular for the configuration exponents 
described in Sect. 2. We can then compare these estimates to the corresponding expected 
theoretical values. 

This section is organized as follows: The case of one infinite river (the original meander 
problem) is discussed in subsection 3.1, while the other geometries, including the case of 
two infinite parallel rivers (connected at infinity), with possibly marked points on each of 
them, and the case of a semi-infinite river (semi-meanders), are discussed in subsection 
3.2. In the three cases above, we allow only crossing vertices of type a (cf. Fig. 1). The 
inclusion of tangency points for the infinite river case (tangent meanders) is discussed in 
subsection 3.3. 
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bridge: 12 3 '". 4 5 

operation: O U O D C* 




Fig. 13: A multi-road meander witli 8 bridges and 2 connected components. 
Upon addition of a new bridge, four different operations (O, C, U and D) are 
allowed, as discussed in the text. We also show the state of that part of the 
meander that is to the left of a given vertical line perpendicular to the river 
(here in position after t — 2 steps). The transfer matrix acts by transferring 
this line from the left to the right. 

3.1. One infinite river (meanders) 

Figure 13 shows a typical meander system with one infinite river and N = 8 bridges. 
The partition function is defined by giving a weight n2 = q to each road. The transfer 
matrix acts by transferring from the left to the right a vertical line intersecting the river 
between two consecutive bridges. A state characterizes that part of the system that is to 
the left of the vertical line by listing the pairwise connectivities amongst the road segments 
as well as the position of the river. 

Whenever the vertical line is transfered one step to the right, a new bridge is added 
and one out of four possible operations (O, C, U and D) can take place, as indicated on 
the figure. Each of these operations connects two consecutive states of the system, before 
and after the addition of the bridge. The operation O opens a new road segment on top 
of the river, thus connecting the two sides of the river. Similarly the operation C closes a 
road on top of the river. This operation comes in two variants, depending on whether the 
closed segment was already connected before the addition of the bridge (C) or not (C*). In 
the former case, the road segment is erased and a non-trivial Boltzmann factor of q must 
be accounted for. In the latter case, the connectivity is transformed so as to connect the 
left-over partners of the two road segments that were eliminated. Finally, a road segment 
immediately below the river can move up (operation U), and a segment just above the 
river can move down (operation D). 
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In order to fully specify the transfer matrix of this problem we need to enumerate 
the possible connectivity states at a given time. The number of such states determines 
the size of the memory needed to store information at time t = n. We will therefore 
compute the number F(n, A^) of connectivity states after addition of the n-th bridge for 
any n = 1,2,..., 2N for the case of 2N bridges. Also, to implement the calculation of 
the meander partition function on a computer, one needs to order these states so that the 
entries of the transfer matrix can be accessed by means of a one-to-one mapping between 
the states and the set of integers 1,2,..., A^states- Here A^states stands for the total number 
of states encountered in the whole transfer process from the left of the first bridge to the 
right of the last one. It is also the dimension of the explored state space and depends 
explicitly on the number 2 AT of bridges added in the whole process. This number will be 
estimated below while the explicit ordering procedure will be presented at the end of this 
subsection. 

(1) (2) (3) (4) (5) 

q 2q +2q 

C".: J IS 1 — 3^-1-3 




^2 
tr^3q+4 



Fig. 14: Transfer matrix enumeration of all multi-road meanders with 6 
bridges. The figure lists the complete set of intermediary states after the 
n-th bridge (n = 1, 2, . . . , 5) along with their respective weight. 

Before we turn to the general case, let us illustrate the state counting for one infinite 
river with 2A^ = 6 bridges. In Fig. 14, we detail the intermediate steps in the calculation 
of the corresponding meander polynomial ms{q). The figure depicts the complete set 
of states generated by the transfer process, when connecting the "empty" states (vacua) 
before bridge 1 and after bridge 6, along with their respective weights. Apart from the 
"empty" states F(0, 6) = F(6,6) = 1, we read F(l,3) = 1, F(2,3) = 4, F(3,3) = 6, 
F(4, 3) = 5, F(5,3) = 1 and A^states = H- The result of the calculation is the meander 
polynomial mz{q) = ^q^ + 12q^ -\- 8q, indicating that with 6 bridges there are respectively 
8, 12 and 5 meanders with 1, 2 and 3 connected components. 



23 



Dimension of the state space: 

Let us now turn to the computation of the numbers F{n, N). To this end we begin 
by relating these numbers to properties of a certain class of restricted Brownian walks. 
Consider a situation where there are pi (resp. P2) road segments above (resp. below) 
the river, and where these segments are pairwise connected in such a way that exactly h 
arches cross the river. It is easy to check that the four operations O, C, D and U described 
above always shift pi, p2 and h by one unit, either +1 or —1, so that, by induction, n, 
Pi, p2 and h have the same parity. Clearly, we also have h < pi and h < p2 as an arch 
crossing the river connects a point above it to a point below it. There are then {pi — h)/2 
(resp. {p2 — h)/2) non-crossing arches that stay above (resp. below) the river. Now, while 
a crossing arch can be generated by means of a single bridge, by using the move O, the 
generation of a non-crossing arch necessitates at least two moves (O followed by either 
U or D). Therefore, twice the number of non-crossing arches plus the number of crossing 
arches cannot supersede the number of bridges added, i.e. pi + p2 — h < n. In terms of /i, 
the arch "height" above the river, this leads to the constraint: 



max(0,pi + P2 — n) < h < min(pi,p2)- (3-1) 



Note that the above condition automatically implies that Pi < n, i = 1, 2. 

The above constraint turns out to be the only one as long as n < A^. For n > N there 
are additional constraints, since we must always be able to annihilate any given state at 
level n in at most 2N — n moves so as to end up with the "empty" state at the right of the 
2A^-th bridge . Since pi and p2 can decrease by at most one at each step, we thus have to 
impose for n > N the two extra conditions: 

Pi<2N-n, p2<2N-n. (3.2) 

It is easy to estimate the number A'^states of accessible states, i.e. those which satisfy 
the above constraints (|3.1|) and (|3.2|) for all n = 1, . . . ,2N . An arch configuration, read 



from the bottom and upwards, can be mapped onto a {pi + p2)-step Brownian walk, 
where the position hi (z = 0, . . . ,pi +P2) of the walk starts from /iq = and is increased 
(resp. decreased) by one unit for each opening (resp. closing) of an arch, as illustrated 
on Figure 15. Evidently these walks are constrained by /i^ > for all i. Decomposing the 
walk into its left part (describing the arches below the river) and its right part (describing 
the arches above the river) connected at height h, we are naturally lead to define /(p, h) as 
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(a) (b) 

Fig. 15: A typical arcli configuration (a) below the river and the correspond- 
ing walk (b). The walk has p = 7 steps and final height h = 3. At each step 
the height coordinate hi on (b) changes by ±1, according to whether an arch 
is opened or closed in (a). 

the number of walks going from height to height h (or equivalently that of walks going 
from height h to height 0) in p steps. 

To evaluate /(p, h) is a standard exercise. First note that without the constraint 
h > this would simply read ((•p_^)/2) must choose the (p — h)/2 ascending steps. 

Each walk violating the constraint at some point hits hi = —1. Let i be the first time this 
happens, and consider reversing each step from i and onwards. The result is a walk going 
from to (—1) — (/i + 1), or equivalently from to /i + 2. Therefore: 

^<-'"K(.-''.)/2)-((.-/>V2-l)- 

In terms of f{p,h), the number F{n,N) of states accessible at level n can then be 
written as 

min(n,N —n) min(n,N —n) min(pi,p2) 

F{n,N)= H E fiPi,h)f{p2,h) 

Pi=0 P2=0 h=ma.x(0,pi+p2-n) 

n min(n,2Af — n) min(n+/i— pi ,2Af — n) 

= E E E f{pi,h)fiP2,h) , 

h=0 pi=h 'P2=h 

where it is implicitly understood that h, pi and p2 in the sums above all have the parity 
of n. 

When iV > n it is possible to rewrite Eq. (|3.4|) in a simpler form, since the second 
constraint in the upper limit of the summations over pi and p2 (namely Pi < N — n) does 
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not come into play. Setting pi + P2 = 2p, we simply have to enumerate all the set of 
p arches, cut by the river at a height h that has the same parity as n and satisfies the 
constraint h > 2p — n. Let us define 

, , ^ f maxfO, 2)0 — n) for n even 

hmin{n,P) = < „ in 11 (3.5) 

max(l, — n) tor n odd. 

In the walk language, we have to count all the walks of length 2p that stay non-negative 
with a marking at a point of height h > /imin and with a well defined parity, that of n. 
The total number of non-negative walks of length 2p going from to and with a marked 
point at position h can be easily calculated to be: 

g{p,h) = f{2p,2h) + f{2p,2h + 2)=(^^)-( (3.6) 



p — hj \p — h — 2 

Summing over the heights h > h^[n{n,p) having the parity of n (which is also that of 
^min(?T',p)) we have 

E 9iP,h)=( y (3.7) 



"mm 

/). = /).jjjjjj(n,p)mod 2 



The complete number of states used in the transfer matrix is now simply obtained by 
summing over the number of arches. Depending on the parity of n, we get: 

^ ( 2p \ 
Fin, > n) = > I I for n even, (3.8) 

^ \Yam{p,n-p)) 

and similarly: 

" / 2p \ 

Fin, N > n) =} { ] for n odd. 

\mm(p — 1, 71 — p)y 



(3.9) 



Note that the above expressions do not depend on A^. This is because we assumed that 
n < A^, in which case the second constraint ( |3.2| ) is ineffective. For large n, the above 
expressions for F(n, N > n) can be evaluated by a simple saddle point approximation. 
Setting p = yn, we get a saddle point at a value y* solution of 



(3y*-l)3 = (l-y*)(2y*)2 , 
26 



(3.10) 



namely with the numerical value y* = 0.611491992 ■ ■ -. This value indicates that the state 
statistics is indeed governed by arch systems with a number p = pi = p2 = p* = y*n of 
arches. The quantity F{n, N > n) is then found to grow asymptotically as a"^, where 



(1 - y*y-y* {3y* - i)^y*-^ 

which, using (|3.10| ), is also the solution of 

a3 = (l + a)2 . 

This yields the numerical value a — 2.147899036 ■ ■ •. 
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21 


7 


134 


161 


8 


19 


1254586 


5733997 


23 


8 


300 


363 


9 


20 


2703503 


14523043 


24 


9 


614 


846 


10 


21 


5777115 


35946838 


25 


10 


1349 


1890 


11 


22 


12438708 


87192966 


26 


11 


2813 


4579 


13 


23 


26613942 


223196395 


28 


12 


6126 


11216 


14 


24 


57268474 


568622062 
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Table 1: Number of intermediate states employed by the transfer matrix algorithm for 
multi-road meanders with 2N bridges. The number of states F{N, N) after addition of 
the A^-th bridge is less than the maximal number of states F{n^i^^, N), which occurs at a 
value rj-max slightly above A^. 

In Table 1, we give some explicit values of the number of states needed when enumer- 
ating meanders with 2A^ bridges. It is seen that the maximum number of states occurs 
slightly after the addition of the middle bridge, n = N . For large A^, one can easily estimate 
the value ?imax of n where this maximum occurs. Indeed, for n > N, the second condition 
( p72|) starts to play a role. As we have seen, for large n, the state statistics without the sec- 
ond constraint is dominated by arch configurations with Pi = P2 = P* = y^n. Therefore, 
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we expect this second constraint ( p.2|) to effectively affect the asymptotic behavior and 
start reducing the number of states whenever p* = 2N — n, i.e. for n = 2N / (1 + y*). As- 
suming that this precisely corresponds to the step having the maximum number of states, 
we have 

?^max = 2A^/(l + y*) = 1.241085907- --A^ . (3.13) 

We can then estimate the asymptotic number of encountered states in the whole process 
to grow like: 

iVstates ~ i^(n™ax, N) ~ a'^— = (2.582603447 (3.14) 
The estimates ( p.l3| ) and ( p.l4| ) agree with the values observed in Table 1. 

Ordering the states: 

According to Eq. (|3.4|) , the full set of F{n, N) states can be ordered if we know how 
to order the /(p, h) restricted Brownian walks defined above. Namely, in that case one 
can order the complete set of states lexicographically after /i, pi, ^2, the value of the 
first Brownian walk [1,2,..., /(pi, h)], and finally the value of the second Brownian walk 
[1, 2, . . . , / (p2, ^)]- By this we mean simply that state A precedes state B ii Ha < hs- 
If Ha = hs this criterion is inconclusive, and one compares the values of pi, so that A 
precedes B if pi^A < Pi,b- In case of further equality one proceeds to compare p2, and so 
on. 

To order the Brownian walks, consider first the example contributing to /(7, 3) shown 
in Fig. 15. The idea is to obtain another formula for /(p, h) which will in turn allow us to 
define a recursive ordering of the walks. We start by focusing on the first (lowermost) arch. 
Either this arch is open (as is the case on the figure), or it closes at some other point before 
p. In the first case the remaining arch configuration is a contribution to f{p — 1, /i — 1), 
and in the second the arches inside the first arch are independent of the arches above its 
termination point. We therefore have 

f{p,h)=f{p-l,h-l)+ Yl f{k-l,0)f(p-k-l,h) , (3.15) 

k — 1 J 3 J 5 J . . . 

with /(O, h) = (5o,/i, and /(p, /i) = for p < 0. This is the required formula. 

A walk contributing to /(p, h) can now be recursively ordered, first by considering 
the termination point of the lowermost arch (which by definition is infinity if that arch is 
open) , then by (recursively) considering the ordering of the smaller arch system inside the 
first arch, and finally by (recursively) considering the ordering of the arch system above the 
first arch (which only exists if that arch is closed). The procedure just described generalizes 



the ordering of the Catalan connectivities Cp/2 = f{p, 0) given in Ref. 
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3. 2. Other geometries 

We now come to the case of the more involved meandric geometries encountered in 
Sect. 2. 

Two infinite rivers: 

It is possible to generalize the multi-road one-river transfer matrix to the case of 
several rivers, provided that the latter can be deformed into a system of parallel lines that 
are only connected among themselves at infinity. For simplicity we consider in the following 
the case of two such infinite rivers. As we have seen, this situation also corresponds to 
the deformation of the figure-eight configuration shown in Fig. 9, provided that a marked 
point is added on each river. 




Fig. 16: The deformation of a two-river configuration by sending all the 
bridges crossing the upper river to the left and all the bridges crossing the 
lower river to the right. The transfer matrix acts by first adding the left 
bridges, then the right ones. 

As before, the crucial point is to have a well-defined transfer direction, which we 
simply take to be parallel to both rivers. For any given configuration with 2A^ bridges we 
write N = Ni +N2, where 2A^i (resp. 2A^2) is the number of times the roads cross the first 
(resp. the second) river. Note that contrary to the bridges on the same river which are 
naturally ordered, bridges on the first and on the second river are not naturally ordered 
with respect to one another. To avoid double counting in the transfer matrix approach, 
we can simply deform the roads as in Fig. 16 so as to send all the bridges of river 1 to the 
left and all those of river 2 to the right. In the transfer process, we thus add first the 2Ni 
bridges crossing the first river and then the 2A^2 bridges crossing the second river. 
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Given a decomposition N = Ni + N2, the calculations therefore proceed exactly as in 
the one-river case, except that when both A^i > and A^2 > not all of the connectivity- 
states described in Sect. 3.1 come into use. Denoting by pi the number of road segments 
above the first river, by p2 the road segments in between rivers, and by ps the segments 
below the second river, we now have the following constraints after the addition of the 
n-th bridge: 

pi < max (min(n, 2Ni — n), 0) , 

P2 < nim{n,2N — n)j (3.16) 

P2 < max (min(n — 2Ni, 2N — n), 0) . 

The total two-river meander polynomial is obtained by summing over the possible 
decompositions. We can furthermore simply address the situation with a marked point on 
either river by weighing each term in the decomposition by a factor of (2A^i + l)(2Ar2 + 1). 

One semi-infinite river (semi-meanders): 

Finally, we have examined systems of semi-meanders, where the roads are allowed to 
wind around the source of a semi-infinite river. Each winding number w can be examined 
separately, by choosing an initial state of w arches nested inside one another, and simply 
applying the four operations discussed in Sect. 3.1. Unfortunately in this situation we 
have not found a simple and full description of the generated states (apart from w = and 
w = 1). It is however still possible to carry out the transfer matrix calculations, by simply 
inserting the generated states in an unordered list. The price to be paid is that to find a 
given state one will have to sequentially search through the entire list, leading to a pitiful 
waste of computation time. This is nevertheless what we have done, and consequentially 
we have had to content ourselves with smaller system sizes. 

3. 3. Tangency points 



OCT 



, ' ' I / ' ^ \ 

C T 

Fig. 17: The six new operations allowed for tangent meanders. 
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We have also enumerated systems of tangent meanders, where apart from the four 
types of moves at each bridge described in subsection 3.1, we aUow for the six additional 
operations depicted in Figure 17. These moves allow to open a road whilst staying above 
or below the river (Oa and Ob), to similarly close a road without crossing the river (Ca 
and Cb), and finally to touch the river whilst staying above it or below it (Ta and Tb). The 
meander polynomial is now defined by assigning a weight x to each of the crossing bridges 
and a weight y to each of the tangent vertices. 

Here again, we have not been able to find an explicit way of enumerating and ordering 
the states obtained by applying this set of ten operations to the vacuum. 



4. Numerical results 

Let us now present our numerical results obtained by use of the algorithms presented 
in Sect. 3, together with an appropriate extrapolation procedure. Subsection 4.1 discusses 
multi-component meanders and introduces the extrapolation method. Subsection 4.2 is 
devoted to the case of two parallel rivers, while subsection 4.3 addresses semi-meanders. 
Finally, we study tangent meanders in subsection 4.4. 

4.1. One infinite river (meanders) 

We have enumerated multi-connected meanders with one infinite river and a fixed 
fugacity q up to N = 24, i.e. 2N = 48 bridges for q = 0, -\/2, and 2. In other words, 
we evaluated the quantities = limg^o(^Ar 5 "^Ar("\/2), mNiV^) and mAr(2) for 
= 1, . . . , 24. From these numbers, we can extract estimates for the large "activity" 
per bridge R{q) and configuration exponent a{q) defined as in Eq. ( |2.11| ) by : 

(4.1) 

The estimates for a{q) can be transformed into estimates for the central charge c{q) through 
the following relation, inverting (|2.17 ) for 7 = 2 — a{q): 



.(,).i_a<i^. (4.2) 

{a{q) - 1) 



Starting from the zero-th order values: 
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we can build better estimates R^^^ and a^'' by a recursive use of a standard convergence 
acceleration procedure: 



R 



ip) 

N 



R 



(p-i) 

N 



p+1 [^N+l " )y^N ~ ) 



P 



p(p-l) I p(p-l) 



2R 



(p-i) 

TV 



(4.4) 



and similarly for a^^^ . The procedure implicitly assumes that the corrections to the asymp- 
totic scaling ( [4.1|) are regular, i.e. integer powers of in which case it guarantees that 
R'^p{q)=R{q) + 0{l/NP+^). 
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Fig. 18: The estimated central charge obtained from the values {q) with 
p = 1, 2, 3 and 4 iterations. The size of the symbols decreases with the number 
of iterations. Each value is represented at an abscissa n corresponding to the 
largest index N — n of mN{q) used in the estimate. We also draw horizontal 
lines at the predicted values of c{q). 

Fig. 18 shows the estimated central charge obtained through ( [1.2| ) using a^p (q) for 
p — 1,2, 3, 4, together with the expected central charge according to (|2.4| ) with ni = and 
77,2 = q. For completeness, we also show the corresponding estimates obtained for q = 1 
with mAr(l) = (cat)^ in terms of the Catalan numbers cn- The corresponding quantitative 
values are displayed in Table 2 below. 
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V2 
V3 
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29+^^ ^ 3.4201328 • 



53+V265 
24 



131 + \/1441 



2.8866175- 
= 2.8160084 



= 2.7675918- 

6 



60 







Ctheor. (?) 


Cnum. (?) 


3.4207 


12.26286 


-4 


-4.003 


3.00000 


16.000000 


-2 


-2.0000 


2.885 


17.52468 


-3/2 


-1.496 


2.812 


18.68970 


-6/5 


-1.18 


2.75 


19.669 


-1 


-0.95 



Table 2: Numerical estimates for the meander configuration exponent a{q), the activity 
R^{q) per pair of bridges, and the central charge c{q). The error is implicitly on the last 
digit. The corresponding theoretical values are also listed. 
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Fig. 19: Numerical estimates (triangles) for the meander configuration ex- 
ponent a{q) for < ? < 2, together with the predicted value ( ^.35| )-( |2.32| ) 
(solid line). We display the result of several iterations of our convergence 
acceleration procedure. The size of the symbols decreases with the number 
of iterations. 



We have also computed the numbers M)^ = meanders with one infinite river, 
2N bridges and a fixed number k of connected components (/c = 1, . . . , A^) for = 1, . . . , 20 
(40 bridges). As an illustration we display the results for = 20 in Table 3 below. 
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k 


^^^20 


1 


64477712119584604 


11 


706958959835806990 


2 


511563350415103008 


12 


235265604762448572 


3 


1901345329566422790 


13 


64713641205591820 


4 


4405839231880790648 


14 


14658557362753320 


5 


7145814923879522986 


15 


2709804590263296 


6 


8632733743310196256 


16 


402058856155712 


7 


8070705247685170684 


17 


46500885666900 


8 


5988061883039308848 


18 


3978168316200 


9 


3587066097601934530 


19 


226760523600 


10 


1755310029771295216 


20 


6564120420 



Table 3: Multi-component meander numbers M^^ with 2N = 40 bridges and k connected 
components of road, k = 1,2, . . . , 20. 

From these values, we can extract the polynomials 'mN{q) = '^^=iM^\'^, hence 
the values R{q) and a{q) for a varying fugacity q. These values are displayed in Fig. 19 
together with the prediction (|2.35| )-( p73^ ). 

All our estimates clearly validate the theoretical predictions with a very good accuracy. 
We note a small discrepancy for values of q close to q = 2. This corresponds to the regime 
which has the worst convergence of our acceleration procedure. This poor convergence 
might be due to either corrections which are not regular, as implicitly assumed, or to the 
vicinity of a transition aX q = 2 where the DPL^ model stops being critical. We anyway 
impute this small discrepancy to our estimation procedure. 

4-2. Two infinite rivers 

Beyond the central charge, we can also test the operator content of the theory. A 
first check concerns the dimension A4 of the operator creating a vertex with four outgoing 
river segments, which, as we already explained, can be measured through the configuration 
exponent 02 -mark for two parallel infinite rivers with a marked point on each river. As 
we mentioned before, in the absence of marking, we expect the two-river configuration 
exponent a2-rivXQ) be simply the same as that of meanders with a single infinite river. 

Fig. 20 and 21 show our results based on enumerations with up to 2A^ = 42 bridges 
for fixed q = 0, V2, \/S and 2. For completeness, we also added estimates for g = 1, as 
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Fig. 20: The estimated configuration exponent for two un-marked rivers, 
obtained witli p = 1, 2, 3 and 4 iterations of the acceleration procedure ( [4.4| ). 
The size of the symbols decreases with the number of iterations. Each value 
is represented at an abscissa n corresponding to the largest index = n of 
m^^^^' (q) used in the estimate. We also draw horizontal lines at the predicted 
values of a2-riv.('?)- 

extracted from the exact expressions: 



2-riv 



for two un-marked rivers, and 

N 



N 

(9 = 1) = XI ^V^N-pCN = CnCn+1 (4.5) 
p=0 



= 1) = p^i2p + l)(2(iV -p) + l)c,CN-,CN = (4^+^ - (^^^^) ) (4.6) 

for two marked rivers. Note that in this latter case of two marked rivers, the leading contri- 
bution m^"™''''^(l) ~ (16)^/A^3/2 is corrected by a sub-leading contribution ~ (16)^/iV^ 
with an exponent displaced by a half-integer. In this case, the corrections to the leading 
term thus involve generically half-integers instead of integers. We expect this phenomenon 
to be generic for all values of q with a double expansion involving both the exponent 
CK2-mark('?) and its "desceudauts" shifted by integers, and a sub-leading correction with 
the exponent a2riv.{Q) ~ 1 descendants. Since a2-riv. — 1 — a2-mark is always close 
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Fig. 21: The estimated configuration exponent for two marked rivers, ob- 
tained witfi p = 2, 3 and 4 iterations of tfie (modified, as explained in tlie text) 
acceleration procedure. Tlie size of tlie symbols decreases with the number 
of iterations. Each value is represented at an abscissa n corresponding to 
the largest index = n of m'^p^™'^^^{q) used in the estimate. We also draw 
horizontal lines at the predicted values of a2-mark(?)- 



to 1/2 for the values of q of interest, we decided to use a modified acceleration procedure 



assuming half-integer corrections to scaling (the index p in (|4.4|) now takes half-integer 



values). The estimates of Table 4 and Fig. 21 have been obtained with this modified 



procedure. 
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Table 4: Numerical estimates for the configuration exponent and the activity R^{q) per 
pair of bridges in the case of two un- marked rivers (upper half) and two marked rivers 
(lower half). The error is implicitly on the last digit. The corresponding theoretical values 
are also listed. 

Here again, the numerical values corroborate our predictions for the various configu- 
ration exponents. Note that we find that the activity per bridge is for these two cases the 
same as for a single infinite river. It should not come as a surprise that the activity -R(?) 
per bridge is independent of the particular meander geometry since it is a bulk quantity 
insensitive to the choice of boundary conditions (river shapes). It is also interpreted as the 
inverse of the convergence radius common to all the generating functions of the numbers 
at hand, also viewed as various correlators within the same theory. 

4.. 3. One semi-infinite river (semi-meanders) 

As a second check of the operator dimensions of the theory, we also considered the 
case of semi-meanders with a single semi-infinite river (ni = 0) and a arbitrary number 
k of connected components of road. Our estimates rely on an exact enumeration of the 
numbers M^^\w) of semi-meanders with bridges, k roads and winding number w for 
l<iV<33, l<A;<iV, 0<M;<iV and w = N mod 2. For illustration, we list 
in Table 5 below respectively the numbers = '^yjM^'^\w) of semi-meanders with 

N = 33 bridges and k roads (irrespectively of the winding numbers) as well as the numbers 
M^^{w) of connected semi- meanders with winding number w = 1, 3, 5, ... , 33. 
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Fig. 22: Numerical estimates for the semi-meander configuration exponent 
a{q) for < g < Qc, together with the predicted value (|2.37| )-( pT32|) . We 
display the result of several iterations of the acceleration procedure. The size 
of the symbols decreases with the number of iterations. 
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21 


1326698396 
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125267070807626 


23 


73827420 


7 


85716694410306 


25 


2638462 


9 


42336073574012 


27 


55052 


11 


15599486790514 
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29 
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13 


4337132101822 


31 


2 


15 


908960663970 


33 





17 


142142103262 







Table 5: Multi-component semi-meander numbers with = 33 bridges and k 

connected components of road, /c 1, 2, . . . , 33, and fixed winding connected semi-meander 
numbers M^\w) with N = 33 and = 1, 3, 5, . . . , 33. 

For this particular geometry, a clear parity effect occurs and it is thus important to 
make numerical evaluations from semi-meander numbers with either even or odd N to 
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avoid large numerical errors. The results presented in Fig. 22 were extracted from even 
values of A^. Table 6 lists quantitative values for the configuration exponent and the 
activity per pair of bridges. Here again, we note that the activity per bridge is identical to 
that of the other geometries. As shown in Fig. 22, the configuration exponent is in very 
good agreement for q less than 1.5. For larger values of a small discrepancy appears, 
again imputable to a poor convergence probably due to the vicinity of a transition point, 
now aX q = Qc where we expect a proliferation of connected components of road for the 
semi-meander geometry. 
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Table 6: Numerical estimates for the configuration exponent and the activity B?{q) per 
pair of bridges in the case of one semi-infinite river (semi-meanders). The error is implicitly 
on the last digit. The corresponding theoretical values are also listed. 

4.4- Tangent meanders 

Finally, we also checked our assertion that the "tangency" vertices of type b (see 
Fig. 1) are irrelevant. We computed the canonical partition function fX]\[{q;x,y) defined 
in (|2.10| ) for the particular case of tangent meanders with a single connected component 
of road {q = 0), a fixed weight x = 1 for "crossing" vertices and with several values of 
the weight y per "tangency" vertex. We enumerated up to = 24 vertices for y = 0.5, 1 
and 2 while y = simply corresponds to the meander configurations of subsection 4.1 for 
which we have results up to 48 vertices. 

Our results are represented in Fig. 23 and the corresponding quantitative values listed 
in Table 7 below. We find a clear evidence that the configuration exponent is independent 
of y, indicating that the universality class of tangent meanders and that of meanders are 
the same, as predicted in subsection 2.3. 
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Fig. 23: The estimated configuration exponent a for tangent meanders with 
a weight y = 0, 0.5, 1 and 2 per tangency vertex. The size of the symbols 
decreases with the number p of iterations (here p = 1, 2, 3, 4). Each value is 
represented at an abscissa n corresponding to the largest index = n for the 
number of vertices (resp. the number of bridge pairs in the case y = 0) used 
in the estimate. We also draw horizontal lines at the predicted value. 
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Table 7: Numerical estimates for the configuration exponent and the activity R{q) per 
bridges in the case of tangent meanders with g = 0, a; = 1 and y = 0, 0.5, 1 and 2. The 
error is implicitly on the last digit. The corresponding theoretical value is also listed. 



5. Discussion and conclusion 

In this paper, we have presented theoretical and numerical evidence for a number of 
exact meandric configuration exponents. Although the sequence of physical arguments 
leading to the exact predictions is by no means a mathematical proof, the numerical 
evidence is compelling. Let us add a few comments on possible generalizations of our 
work. 



41 




Fig. 24: The range of validity of ttie exact prediction for the meander 
configuration exponent in the general case of a weight ni per river and n2 
per road. In addition to the criticality constraint that < ti^ < 2, we have 
represented the c(ni,n2) = 1 curve (joining the points (tii = 2,n2 = 1) and 
(ni = 1,77-2 = 2), and passing by (ni = \/2, n2 = a/S)), beyond which the 
meanders are in a branched polymer phase. We have also represented by a 
dashed line the location of the winding transition joining the point (0, Qc) to 
the point (1, 2). 

5.1. Ranges of validity and transitions 

The physical arguments presented above allow to go much farther than just the case 
of one or two rivers. Indeed, we have seen that meanders are a particular case of the 
GFPL^(7ii, 772) model at rj-i = 77,2 = with y = 0. Let us define multi-river and multi-road 
meander polynomials of order 71, 771^7,(77.1, 77-2), by the following expansion of the partition 
function -Z'gfpl ( |2.5| ) for ?/ = 

^GFPL(ni, 712; x) = -— m„(?7i, 772) . (5.1) 
n>l 

Then we can as well predict the following large 77 asymptotic behavior 

R{n\ 772)^*^ 

"^n(ni,?72) ~ C(?7i,?72) — ]l l , , a(ni , 772) = 2 - 7(c(?7i , 7^2) ) , (5.2) 

where c(ni,n2) is the central charge of the DPL^(ni,n2) model (p.4|), and 7(c) is as in 
( p.lTp . The general prediction ( [5.2| ) was actually proved in the case 72,1 = 1 and arbitrary 
< 772 < 2 in [0, by solving a particular matrix model. 
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Note however that the range of vahdity of ( |2.17| ) imposes that c(ni,n2) < 1. This 
gives the total range of vahdity of Fig. 24. Outside of this range, we must consider the two 
foUowing cases: 

(i) One of the Ui > 2: then the DPL^ model is no longer critical, and we have no field 
theoretical prediction as to the value of the configuration exponent. It has been shown 
however, using matrix model techniques, that in the resembling case of the so-called 
0(n) model, the critical exponent takes a constant value a = 3/2 independently of 
n > 2 0. 

(ii) Both ni,n2 < 2 but c(ni,n2) > 1: one encounters the well-known "c=l barrier" 



phenomenon in two-dimensional quantum gravity |]31[]. The corresponding theories 
are dominated by configurations of surfaces with long fingers (branched polymer phase 
of quantum gravity). Remarkably, it was shown that 7 = 1/2 throughout this phase, 
leading also to a constant exponent a = 3/2. 
Hence we expect (for different reasons) that a = 3/2 identically outside of the range of 
Fig. 24, for ni,n2 > 0. On the ordinate axis of Fig. 24, where ni = 0, we recover for 
meanders the announced transition at 77-2 = Q = 2, beyond which the DPL^ model is no 
longer critical; numerical results however get poorer as we approach this point and we 
could not confidently analyze it numerically so far. In the case of semi- meanders, we have 
predicted an earlier "winding" transition point at n2 = Qc = 2cos7r(v^97 — l)/48). More 
generally, the configuration exponent for multi-component and multi-river semi-meanders 
(one half-line plus an arbitrary number of river loops) is predicted to be 

a{ni, 77.2) = a(ni, n2) - 1 + 2Ai(ni, 77-2) 

1 (5-3) 
= 1 + ^ ( ^25 - c(ni, 772) + - c(ni, 772)) ^6(1 - ei) - 4c(n2) , 

where 0(772) = 1 — 6e|/ (1 — 62) and rti = 2 cos nci. The range of validity of ( |5.3[ ) is smaller 
than the domain of Fig. 24, as it is delimited by the curve 3(1 — ei)/2 -(- 6e|/ (1 — 62) = 1, 
with Ui = 2cos7rei (it is represented in dashed line on Fig. 24). The latter corresponds 
to a winding transition as explained above, where the number of pieces of road winding 
around the source of the semi-infinite river becomes relevant. Above that curve, we expect 
the semi-meander configuration exponent to identically vanish. 
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5.2. Extensions 

Extending the generalization of the previous subsection, we could consider both com- 
plicated river geometries and keeping ni,n2 finite. This requires some care when dealing 
with the electromagnetic operators creating river vertices, as several geometries might 
correspond to the same correlators. Also extra care should be exercised when imposing 
that the meandric objects be connected: in the multi-river and multi-road case, a meander 
must be globally connected, but can have disconnected rivers or roads. We may also de- 
fine higher-genus meandric numbers ^ which can also be enumerated by transfer matrix 
techniques, for instance on a torus. 

As discussed before, we may also consider the coupling of the FPL^ model to Eulerian 
gravity, namely by summing over genus zero tetravalent graphs with only faces of even 
valency. Then all the above predictions are expected to still hold, except that we must 
take the formula ( p.3| ) for the central charge, that remains one unit above that of the DPL^ 
model as we have ensured the vertex-bicolorability of the graphs we sum over [^. Note 
that, in this Eulerian case, the b- vertex of Fig. 1 is now relevant and we must take y 7^ to 
get the above shifted central charge. The range of validity of Fig. 24 is reduced in this case 
to the zone of the square [0, 2] x [0, 2] delimited by the curve c = 1 joining the point (2, q*) 
to the point (g*, 2), with q* = 2 cos(7r(v^- l)/6) = 0.410135 • ■ -. One should also be able 
to check numerically the predictions for Eulerian gravity, using a suitable modification of 
our transfer matrix technique, in order to incorporate the vertex-bicolorability. 

Another direction of generalization would be to consider more loop colors. For in- 
stance, one can define folding problems of two-dimensional lattices onto a (i-dimensional 
target space, that in turn correspond to fully-packed loop models with a given number 
of loop colors, related to d. Attaching a different fugacity for loops of different colors, 
we may generate different universality classes, described by different conformal theories. 
Upon coupling to two-dimensional quantum gravity, these models would correspond to 
multi-colored meanders of some kind. 

A final direction of generalization would be to add "matter" to the meandric graphs. 
Indeed, imagine we would like to consider a more involved model for a compactly folded 
polymer (protein), by attaching to its nodes a spin variable with inter and intra-chain 
Ising-like interactions. Then a simple way of describing it would be by first defining this 
matter spin model on the configurations of the square lattice fully-packed loop model, and 
then switching on gravity. If the matter model is still conformal, this should immediately 
lead to new configuration exponents for meanders in the presence of matter. 
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5.3. Algebraic exponents 

It is both interesting and sad to notice that the only meandric numbers that can be 
calculated exactly with reasonably simple combinatorial formulas are all in one way or 
another related to Catalan numbers, and the corresponding exponents are always integers 
or half-integers (take for instance the case ni = 0, 77-2 = 1, then c = —2, a = 3, a = 3/2). 

For general rational values of the e^, we predict however that the exponents are alge- 
braic numbers, roots of some quadratic equations with integer coefficients. Remarkably, 
such exponents are not commonplace in physics, but have emerged in some recent works on 
random walks . If we could devise some relation between the meander problems and 

this other type of problems, we would probably be able to get a mathematically rigorous 
proof of our predictions. 

Acknowledgements: We thank O. Golinelli for useful discussions. 
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